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TRANSFORMATIONS OF GAUSSIAN LIGHT BEAMS 
CAUSED BY REFLECTION IN FEL RESONATORS 


I. INTRODUCTION 

Free Electron Lasers (FELs) operating as oscillators* -2 require the 

8 9 

trapping of light pulses between systems of mirrors (resonators). * These 
pulses are repeatedly amplified via synchronous interaction with electron 
pulses passing through the wiggler. The radiation produced by the 
stimulated emission is confined within a narrow cone along the beam axis. 
Therefore, the vector potential can be represented as a superposition of 

i kz 

Gaussian modes. Those are the free space eigenmodes A mn (r) = e mn A mn^ r ^ e 
where e is the polarization vector, of the paraxial equation,*® 

V x 2 A - 2ik §f = 0. (1) 

Equation (1) is the k^ « k = w/c limit of the exact wave equation. The 
simplest oscillator configuration is that of an open resonator with two 
opposed identical mirrors. The vacuum eigenmodes for this arrangement are 
also expressed in terms of the paraxial eigenmodes. Their detailed 
structure can be described in terms of either Gaussian-Hermite functions in 
rectangular coordinates, or Laguerre functions in polar coordinates. In 
both representations all the eigenmodes with given wave number k are 

1/2 

characterized by two independent parameters: the waist w = (2b/k) and 

2 2 

the curvature of the wave front 1/R = z/(z + b ), where z is the distance 
from the waist position and b is the Rayleigh length (Fig. 1). 

The electron beam is an optically active medium that alters the 
characteristic parameters of the radiation after each passage. During the 
build-up period the modal content and the structure of the light pulses 
inside the oscillator will change. A numerical method has been developed 
recently optimizing the representation for the amplified radiation. In the 
11 12 

source dependent expansion ’ the waist size and the curvature of the 
elected modal basis is tailored according to the driving source term. That 
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minimizes the number of modes required to describe the light beam. In 
general, the curvature and waist size of these modes does not match the 
curvature and waist of the vacuum eigenmodes for the resonator. Therefore, 
the transfer matrix for a given mirror must be known for arbitrary incoming 
modes. This need stems from computational as well as physical reasons. 

The knowledge of the cavity reflection matrix R, together with the gain 
matrix G through the wiggler, is necessary in determining the potential for 
steady state operation. 

The study of the reflection matrix must include the effects of 
deflecting the light beam, in addition to finite mirror size and curvature 
mismatches. During high power operation, grazing mirror incidence may be 
necessary to avoid exceeding the dielectric breakdown limit for the 
reflecting surface. Also, in case of a high per-pass gain with optical 
guiding, the spot size for the reflected radiation could be much larger 
than the incoming. In two mirror resonators, the reflected radiation could 
then damage the wiggler. Therefore, ring resonators, including three or 
more mirrors, must be employed for the deflection and recirculation of the 
light pulses. 

We are interested in cases when the reflected radiation remains 
focused along some direction z q making an angle ♦ with the incoming z^. 

Then the reflected vector potential will also be expandable in free space 
eigenmodes Ap^(r Q ) of the paraxial equation in the new direction. The 
mirror surface generating focused reflection in the desired direction can 
not be arbitrary but must be appropriately defined. The angle of 
deflection ♦ will enter the equation defining the mirror surface. The 
other surface parameter, namely the curvature 1/R m , is a free parameter. 

It determines the curvature 1 /R q for the outgoing modes given the curvature 
1/Rj of the incoming modes. In case of reflection by an arbitrary surface, 

2 







the scattered radiation cannot, in general, be covered by the paraxial 
modes that do not form a complete set in three dimensions. 

A single incident mode A ran ( r j) will* in general, be partially 
reflected into different modes A pq( r Q ) where (m,n)*(p,q). This is caused 
by the deflection of the light beam, the finite size of the mirror and the 
curvature mismatches. Reflection into other modes will affect the 
interaction between the electron beam and the radiation in a number of 
ways. First, the rms radius of the light beam will change, affecting the 
matching beam condition. Second, the light pulse will spread axially 
because of dispersion among different modes, since the phase velocity 
depends on the modal number (m,n). Third, different phase shifts among the 
various modes during reflection may render these modes out of phase after a 
number of bouncings off the resonator. For the above reasons the fraction 
of radiation scattered into other modes will contribute to the losses in 
FEL oscillators. 

The method for obtaining the reflection matrix is outlined in Sec. II. 
The definition of the appropriate mirror surface is given in Sec. III. In 
Sec. IV the integral expressions for the matrix elements are derived. An 
analytic expansion in powers of a small parameter (of the order of the 
diffraction angle) is given in the same section. Some limiting cases are 
examined in Sec. V. In Sec. VI the reflection of the fundamantal mode 
(0,0) is studied in detail. Section VII deals with cross-coupling effects 
among the vector components of the radiation. 


3 





II. OUTLINE OF THE METHOD 


The free space eigenmodes A mn (r) of the paraxial wave equation have 
the general form 




( 2 ) 


The first exponential in (2) contains the rapidly varying phase on the 
wavelength scale X = 2rt/k. The wavefronts are spherical with radius of 
2 2 

curvature R(z) given by 1/R(z) = z/(z + b ). The spot size of the 
radiation envelope is W(z) = w(l + z^/b^)*^, where w = (2b/k)*^ is the 
waist, and the distance z is measured from the position of the waist. The 
amplitude squared of the mode drops by 1/2 over a distance equal to the 
Rayleigh length b (also known as confocal parameter). Most of the 

radiation is confined within a cone parametrized by the diffraction angle 
1/2 

9j = W/z = (X/bn) . The structure of the amplitude profile u mn ( r >W) 

depends on the elected coordinate system. u mn (r;V) contains the slow 

spatial variation equivalent to a small wave number perpendicular to the z- 

direction. Higher modes correspond to an increasing effective k^, 

producing the slow phase shift expressed by the term exp [iS mn (z)]. For a 

given k, the mode is completely defined by the two independent parameters R 

and w (or any combination of two out of the four quantities R, w, z and b). 

The geometry of the reflection is illustrated in Fig. 2. The 

subscripts i and o denote the coordinate system used for expressing 

incoming and outgoing modes, r.. is defined with the z^ axis along the 

direction of incidence and r has the z axis along the direction of 
o o 

reflection. The origins are displaced from the mirror center by lj and 1 Q 
respectively, where 1^ is the distance of the minimum waist w,, for the 




incoming radiation and 1 Q is the distance of the minimum waist w q for the 
reflected mode. A third coordinate system r with the origin at the mirror 
center and z aligned with the mirror axis will be useful in the 
computations. Underlined quantities rj, r Q and r g stand for the mirror 
surface coordinates in each reference frame. The transformations among the 
various frames are defined by 


z sin 
s 2 


y i 


(3a) y 0 = y s , 


(3b) 


4> 

cos £ + : 


;in | i 


Me consider incoming radiation of given curvature and of arbitrary 
amplitude profile A^Tj), consisting of various modes (m,n) with the same 
Rj(z). If both incident and reflected radiation are expanded into 
eigenmodes, 


A 1 (r.) 


'W v i 

e £ c mi 


it (r ) 

.O/ x o o r» 

A (r ) = e L ' 

p.q 


(4b) 


where 


2R.(z) 


5 


*.(r) = k z 


(4c) 





the relation among the incident and reflected expansion coefficients c' 


is formulated by 

o » i 


i 

mn’ 


(5a) 


or 


o 

pq 


E 

m,n 


R mn 

pq 


(5b) 


where R mn are the elements of the reflection matrix R. 

pq 

We examine the case when the mirror dimensions p are much larger than 

the wavelength X, X « p (otherwise diffraction rather than reflection 

would prevail). We also assume that the angle C, subtended by the mirror 

C = p/R , where R parametrizes the radius of curvature, is small, of the 
m m 

order of the diffraction angle 0^, C ~ 0^ ~ t. The v-th component of the 
reflected vector potential at distance |r Q - r Q | » X from the mirror 
surface S is then given by 


(v) v 


. ik [[ £ 

2n £ 


e ik IV Hoi 
k - Lj 


(r ) (n • Sr). 


In Eq. (6) n*Sr is the obliqueness factor where Sr = (r Q -r o )/|r Q -r o | and n 

is the normal unit vector to the reflecting surface. The surface element 

ds is given by ds = 5[z -f(x ,y )ldx dy dz where z = f(x ,y ) is the 
° * o o ■'o * 0^0 0 O 0 ^0 

surface equation. Equation (6) is the convolution of a source term A s (r ) 
at the mirror surface with the propagator exp(ik|r Q -r o |)/|r Q -r o |, i.e., a 
superposition of spherical waves originating at S. The source term A s (r Q ) 
is specified from the incoming vector potential A*(r.) through the boundary 
conditions and the coordinate transformations (3). We will assume a 
perfectly conducting surface, where the incident and reflected fields are 
related by 

6 







A S = -A 1 + 2 (n • A 1 ) n, 


(7a) 


and n is the normal unit vector to the reflecting surface. The second term 
in (7a) introduces a coupling among different vector components, caused by 
the mirror curvature. This cross-coupling is small and disappears in the 
plane mirror limit, 


s - - A 1 
(\>) " A (%>)’ 

. i 


(7b) 


where A and A are expressed in the incoming and outgoing coordinate 
systems respectively. Because of the linear superposition principle, Eq. 
(6), the cross-coupling contribution can be added separately, and will be 
deferred until Sec. VII. In the next three sections we will treat the 
reflected vector components as independent scalars, according to (7b), that 
corresponds to a phase shift by it during reflection. Most of the 
computations will be performed on the mirror surface. To simplify the 
notation from now on, we drop the bar (_) under the mirror coordinates r. 
Subscripted quantities such as r^, r Q , r g will signify the mirror 
coordinates in each reference frame. Unsubscripted coordinates will denote 
the observation point in the reflected radiation frame of reference. 

We seek cases when the reflected radiation propagates along z q , 
contained within a cross section of dimensions x,y « z-z q . The expansion 
|r-r Q | = (z-z q ) (1 + [(x-x q ) 2 + (y-y Q ) 2 J/2(z-z q ) 2 } replaces the full 
propagator inside (6) with the paraxial propagator in that direction, 


'(r) = JJds A(r Q ) (n • 6r) U k (r,r Q ), 


(8a) 






where 


ik 

2n 


(x-x ) + (y-y ) 


(8b) 


Expression (8) is the approximation of the exact solution (6) to order 
[(x-x q ) 2 + (y-y o ) 2 ]/2(z-z o )^ ~ e 2 . It is valid provided the surface S 
produces focused reflection along the desired direction. Otherwise the 
paraxial limit will fail to encompass all the radiation contained in the 
original expression (6). The geometry of the mirror that is compatible 
with the above approximation will be obtained during the computation of the 
reflection matrix. 

It is known that the profile of a given eigenmode A mn ( x 0 *y o > z 0 ) at z 0 
is generated by the propagator U^r,^) acting on the mode A mn (x,y,0) at z 
= 0. The inverse propagator U k (r,r o ) therefore reproduces A mn (x,y,0) from 

wvW’ 




) = u (x,y,0). (9) 


This suggests multiplying and dividing the integrand inside (8a) by 
2 2 1/2 

exp[i*(r Q )l / (1 + z q /b Q J , recasting (8a) in the form, 

o rr iA(r o ) - i * 0 (r o ) 

A°(r) = JJds e 0 S(r Q ) e U_ k (r,r Q ), (10) 

where the source term S(r Q ) is, 

[ z 2 (r )1 
1 + ° b 2° 1 
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(11) 







and the phase term A(r ) = $.[r.(r )] + * (r ) is given by, 


z i< 


x. (r ) + y. (r ) 


2 2 « 

> + y o 

V r o> .‘ 


The phase A(r Q ) depends on the angle <j> through the coordinate 
transformations Eqs. (3). 

The term exp[id(r Q )J is varying rapidly, on the scale of the 
wavelength X. Therefore, its convolution with the slowly varying source 
term over an arbitrary surface will be vanishingly small. In general, this 
corresponds to radiation scattering where only a small fraction of the 
incoming radiation is reflected along the considered direction <J>. The 
integral (10) will be finite only when it is possible to satisfy the 

condition A(r ) = constant over some surface S. If, in addition, S is much 

o 

larger than X, expression (10) will be finite only within a narrow angle 6$ 
around <t>. This guarantees that the reflected radiation remains focused 
along that direction. Therefore, a condition that the exact reflected 
radiation (6) be fully covered by the paraxial limit (10) is that 


A(r Q ) = constant, 


(13) 


along the surface S. Accordingly, the optical path is the same along the 
rays connecting an incoming wave front with its mirror image (reflected) 
wave front. 

Requirement (13) defines the appropriate mirror surface z o =f Q (x o ,y o ;♦) 
for reflection in the elected direction. Assuming that f Q is found, we may 
express z q in terms of x q , y Q and use the constancy of A(r Q ) over S, 
reducing (10) into 


dx dy o(x ,y ) e 


9 


A°(r) 


U -k (r ’ r o>- 


(14) 




a(x o ,y o ) s s t WVVV* is full y expanded in terms of u mn ( x 0 *y o ) that 
form a complete set in tvo dimensions, 


o(x.y o ) = E Rml 


(15) 


The expansion coefficients R mn for Gaussian incoming radiation of 
arbitrary profile a(x Q ,y Q ) are given by 

R ” - JR, d ’’o «<W u mi/ x o’^o’ R o^ JJ dx o dR o u mn^ x o ,R o’ R o^ ‘ < 16 > 

The radiation spot size W q at the location of the mirror center is a free 
parameter, yet to be specified. Each choice of W Q generates an equivalent 
representation for «r(x o ,y Q ). 

Upon substituting expansion (15) inside the integral (14) and using 
the property (9) for the inverse propagator U_ k , the reflected vector 
potential assumes the final form 


A°(x,y,0) = Z R mn 
m,n 


(17) 


? ? i n 1/? 

where W q (z) = w q (1 + z z /b Q ) , w q = (2b Q /k) . Expression (17) is a 

complete decomposition of the reflected radiation into paraxial eigenmodes 
for incident radiation of arbitrary profile. Therefore, condition (13) 
that defines the mirror surface is sufficient for the full reflection of 
paraxial (Gaussian) incoming light beams into paraxial beams only. The 
fraction of the electromagnetic flux incident on the mirror is conserved 
after reflection. If, on the other hand, (13) is seriously violated, the 
paraxial modes are inadequate to include all reflected radiation, and the 
incident flux is not conserved by expressions similar to (17). 
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III. MIRROR SURFACE. 


To obtain the equation for S ve express all quantities inside (12) in 
the mirror coordinate frame applying the transformations (3a) and (3b). 
Using the scaling x s /R m - y s / R m ~ e « 1, Z s / R m ~ e 2 we obtain from (13) 

[x g 2 cos 2 | + y s 2 ], (18a) 


Equation (18a) is the analytic expression for a paraboloid surface. 

R m parametrizes the mirror curvature, being positive or negative for a 
convex or concave mirror respectively. The surface is reflection symmetri< 
with (zx) g and (zy) g as the symmetry planes; there is no rotational 
symmetry around z g . Surface (18a) can also be approximated, to second 
order in (x /R ) 2 , (y /R ) 2 by hyperboloids or ellipsoids defined by 


.2 ± v 2 

’ 2 y s 

= R ra 2 cos 2 |, 

(19a) 

.2 4 + y 2 

’ 2 + y s 

= R m 2 cos 2 |. 

(19b) 


All the surfaces become spherical in the limit of perpendicular incidence $ 
= 0, and plane mirrors when R m ■* ®. Using the definition of the curvature 
for the paraxial modes, Eq. (2), and the fact that R » b in cases of 
interest, ve obtain from (18b) 

1 2 1 


Relation (20) defines the curvature of the reflected modes from the 
incoming mode curvature and the curvature of the mirror. 





Equations (18)-(20) imply that 


A(r s ) a = const. + 0 [ kp [jp] 1, (21) 

m 

where p parametrizes the mirror size. A more complicated surface equation 
(higher than quadratic in x, y, z) is required to improve the constancy to 
a higher order. In the next section the reflection matrix will be computed 
by expansion in powers of w 0 /,R m - p/ R m - Si nce kp » 1, the approximation 
A(r g ) = constant is satisfactory for a first order expansion as long as 
p/R^ - 1/kp. In case that p/R m > 1/kp, A(x g ,y s ) is a slowly varying 
function over S. Large mirrors require the inclusion of the phase slippage 
term exp [iA(x s> y s )J next to the source term o(x s ,y g ) in Eq. (16). 

The unit vector n normal to the mirror surface is given by 


where 


£ (x g , y g , z g ) is given 


by Eq. 


(18a). 
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IV. COMPUTATION OF THE REFLECTION MATRIX 
According to the definition (5b), the R™|j element of the reflection 
matrix R is obtained from the source term ff pq( x 0 'y o ) inside (14) generated 
by a single incident eigenmode A p( ^[r^. (r Q > ]. The integration is performed 
in the mirror-aligned coordinates, taking advantage of the existing 
symmetries. The coordinates r^ and r Q , defining the incoming and outgoing 
wave functions, become explicit functions of x g , y g through the 
transformations (3). The surface equation (14a) is used to express z in 


2 2 4 2 

_ cos f + y„ 


is defined by the intersection of the infinite surface (18a) with the plane 


reflection matrix elements take the form 


u mn ( x _ * y c ) u n n( X c * ) 




i6(x,y) - i6_(x ,y ) . , x„ , x„ sin £.1 

p<! s s "" s s cos f (l - f sin * - 5 2 2 ) ,(23) 
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2 2 

Expression (23) is correct to order p /R^ . 

It will be seen that R, as given by (23), depends on four parameters 

R s R(4>, «, y; £). (25) 

$ is the reflection angle shown in Fig. 2. a is the ratio of the incoming 

to the outgoing spot size at the mirror, a = W.(1 j)/V o ( 1 q ). y = p/W Q 

parametrizes the mirror size compared to the radiation spot size. £ = 

W /R scales as the diffraction angle 0, = V /I multiplied by the 
om d o o 

curvature mismatch R /R between the mirror and the radiation wavefronts, 
o m 

The spot size W q enters as a free parameter because only the curvature 1 /R q 
for the reflected modes is specified by the mirror geometry. Since many 
combinations of W Q and 1 Q apply to a given curvature according to paragraph 
Eq. (2), an additional selection rule for W q is needed. Note that W q does 
not have to match VT. This is obvious in cases when the mirror size p is 
smaller that W.. Each value of W q defines a complete set of modes for the 
reflected radiation and an equivalent representation for R. 

Parameters <j>, a, and y can be arbitrary. In most cases of interest, 
however, £ is small, £ « 1, of the same order as the diffraction angle 9^. 
The analytic computation of the matrix elements is carried out by expanding 
the integral (23) in powers of £, 

R = R(0) + Z R(1) + l 2 R(2). (26) 

Each representation of R is tied to the choice of the basis functions 
u mn (r). The eigenmodes u mn (r) are specified according to the coordinate 
geometry. In the next subsections we derive R in Gaussian-Hermite and 
Gaussian-Laguerre representations. For simplicity, it is assumed that 
A(x g ,y s ) in Eq. (23) is constant, i.e., kp ( p/R m ) << 1. 
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(a). Gaussian-Hermite representation 

In rectangular coordinates (x,y,z) the functions u mn (x,z;V) are 
given by 


»„(#9 


-W 

„2 


(27a) 


where H , H are the Hermite polynomials and a is a normalization factor, 
m n mn 

setting the total electromagnetic flux carried by the mode equal to unity, 


/2 ( 9 m- 

■» t” 2 


(27b) 


The corresponding slow phase factor & mn ( z ) in Eq. (2) is 


6 mn (z) = (m ♦ n + 1) tan" [— j. 


Substituting inside (23), expanding in £ and performing the integrations, 
Eqs. (23)-(26) yield 


R mn (0 ) = c mn g pq jmi 


R mn (l) = C mn e pq 


where C is a normalization factor 


(28b) 


! /m" 10 + i N 1 " 11 ) , 

\ pq pqj 

actor 

w , .-1/2 

C™ = w: (2 m+n+P+q m!n!p!q!j , (28c) 

the phase is expressed by 

= (p+q+Dtan -1 ^) - (m+n+Dtan" 1 (j^) + k(l i+ l o ). (28d) 
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JdX JdY H p (aX) H q (Y) HJX) H n (Y) , 


s s 
X Y 


. JdX JdY R p («X) H q( Y) H„(X) H„(Y) , 


<* +1 /v 2 . v 2s 


« +1 , v 2 . v 2 x 


f I x * ^ x(x2 * ?2 > - (« b^tSt - irS") r < x2 * v2 >} 


(29b) 


N mn = (dX fdY— H (otX) H (Y) H (X) H (Y) < 

pq j j /£ p q m n 


directly. The matrix elements are computed, to first order in £, in 
Appendix A for large size mirror and a = 1. 


(b). Gaussian-Laguerre representation 

2 2 1/2 

In cylindrical coordinates (r,6,z) where tan© = x/y, r = (x +y ) , 

u^(r,9;W) take the form 


u^<r,9;y) - aP 


fcosp©' 

vsinp©, 



(30a) 


where +p(-p) signifies cosine (sine) poloidal dependence, a^ is given by 



(30b) 
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and the L** are the Laguerre polynomials. The corresponding slow phase 
i£(z) in Eq. (2) is, 

$P(z) = (2m + p + 1) tan" 1 (30c) 


The transformations among polar coordinates representing the various 
reference frames are 

r^ a r g [l - sin 2 0 g sin 2 | - 2 ^ sin 0 g sin | cos 

r o = r s I 1 " sin 2 0 g sin 2 | + 2 ^ sin 0 s sin f cos f] > (31a) 


tan 0. = cos ? tan 0 - — 


. z sin I 

2 tan 6 s " cos 0 


* z s sin 2 

tan 0 Q = cos f tan 0 g + - ^~ r . 


(31b) 


The mirror surface (18a) is expressed in polar coordinates as 

r 2 sin 2 0 cos 2 I + cos 2 0 
_ s _s_2_s 


Applying similar computational procedure as in the previous subsection we 
obtain 


- C” e i * P \o S t 


; { « fin* B P "(X)[v^(X).i V*. ] 


17 




_pq _ 1 _ r mini a 2 1 1/2 
mn ' 2n Um+p)! (n+q) ! J 


and the phase is expressed by 

= ( 2 m+p+l) tan' 1 ^) - ( 2 n+q+l) tan -1 ^) + k(l.+l o ) • ( 32 c) 


D -,x, .[“a. C "A < V]”?J5. < V] 

o J m 1 - sin j sin ^ e m 

r 2n sin 0 cos p[e. (0 )] cos q fe 

E u> ’I 


_ sin 9 m cos P L 0 i <e s ) J COS q LWj 

(1 - sin 2 * sin 2 0 m ) 3 ' 2 

E 9. _ 2L±1 v 

■ 2 x^ 2 x 2 i. p r« 2 x^ i. q rx^ p 2 


U PM (X) = (a X ) 4 if L H (a X) I/VX) 




( „ 

2 i 


fI P (a 2 X) 

L q (X)l 'l 


1 

* /2 

{*?-* 

a +1 
" 2 

X + 





E 

a±i 



2 1 
a +1 



(«V 

X 2 


X) L q (X) e 

" 2 x 



„ E 

a±i 



2 . 

«_li x 

“ S«) 

1 

✓2 

/ 2 V .2 v 
(a X) X 

2 L P 
m 

(« 2 X) 

L S<*> e 

2 

In obtaining 

(33a) 

- (33e), 

X was 

def i 

ned by X = 

ll-sin 2 (+/2)! 

thus, accordi 

ng to 

(22) and 

1 (30), 

the 

boundary 

X is X = 
s s 

lowest order 

terms 

for the 

first 

few 

elements 

are given in i 

arbitrary deflectii 

on angle 

$ and 

a = 1 




;in 2 0 ]r 2 / 2 W o 2 ; 
- 2 /ff 2 
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V. LIMITING CASES 


When the mirror radius tends to infinity (1/R m ■* 0), or in cases of 
vertical incidence on the mirror ($ = 0), the higher order corrections in 
the reflection matrix R disappear, 

R = R(0) (34) 

in both representations. The nondiagonal elements in R stem from the 
finite mirror size only. If, in addition, the mirror size is very 
large, u » 1, it is appropriate to take W Q = as best representation 
for the reflected radiation. The a = 1 limit yields 



Thus, in case of large curved mirror and vertical incidence, or large 
plane mirror and arbitrary incidence, the reflection matrix is the 
identity matrix. 

The case « = 1 is of special interest for arbitrary angle of 
deflection ♦ and mirror curvature 1/R, as it will be explained in the 
next section. For finite mirror size p > W q , ( u > 1), there exists 
zeroth order non-diagonal terms inside R(0). Since R(0) is independent 
of the angle of deflection $, the finite mirror size yields the dominant 
contribution to the reflection into modes different than the incoming. 
The effects of the deflection of the light beam enter to first order in 
£., R(l), or higher. In the Hermite representation the elements Rp^(O) 
couple mode combinations with m + p = even, n + q = even. The elements 
with either m + p or n + q odd vanish because of the even/odd symmetry 
of the Hermite functions. 

As the mirror size becomes very large and the limits of integration 
in (23) are extended, the orthogonality among the various modes u^ (r s ) 
becomes effective. The off-diagonal terms in R(0) become comparable to 
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2 

the first order corrections roughly when 1/p ~ £ ~ 0^. At the limit u 

-» ® all the nondiagonal elements of R are reduced to order £ or higher, 

R pq = * R pq< X > + °<^>» ■ 96 P» n * q, (36a) 

and the only matrix elements of zeroth order in £ are the diagonal 

»In ' C<°> * < 36l » 

in both Hermite and Laguerre represenations. The lowest correction in the 
2 

diagonal elements is of second order £ , while the first order contribution 
disappears. This is consistent with flux conservation during reflection in 
case of large mirror. 

In obtaining Eqs. (28) and (32) it was assumed that A(x g ,y g ) is 
constant over S. According to (21) the variation of A is parametrized by 
£.* = (kV^/R m ) l. When (kW.^/R m ) >1, becomes of order £ and the 
effects of the slow phase slippage must be retained in (23). This effect, 
known as spherical aberration, causes additional corrections R (1), of 
order £, , 

R = R(0) + E R( 1) ♦ £. R (f) + . 

Spherical aberration does not disappear at the limit of large mirror 

size, as opposed to the effects discussed so far. In fact, when £,* > £, 
it places a lower limit on the off-diagonal terms in the reflection 
matrix, 


wm pq * 

Perfect reflection, requiring E. = 0, is possible only for plane mirror (R m 
-> <°) of large size. 

The superposition principle can be used to describe reflection from 
more complex mirror surfaces. In case of a mirror with a hole the surface 
integral (14) over S m is expressed as /<. = where is defined by 

rror exterior boundary and S 2 is the surface of the hole. The total 


the mi: 






reflection matrix R is given by R = R(Sj) - R(S 2 ), the difference in the 
reflection matrices associated vith mirrors Sj and S 2 respectively. The 
transmission matrix T through a screen vith an aperture of area S is given 
by T = - R, R being the reflection matrix for a mirror matching the 
aperture S. The transmission matrix for radiation diffracted behind a 
finite size mirror is given by T'= 1 - e 1 R where 1 is the identity 
matrix. 
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VI. REFLECTION OF THE LOWEST ORDER MODE 


The computation of all the truncated integrals for finite mirror 

surface is nontrivial. Most applications, however, involve the (0,0) 

lowest order mode as the dominant mode in both incoming and reflected 

radiation. The strategy here is to compute the element Rqq of the 

reflection matrix first. Then the waist for the reflected modes W can be 

o 

00 

selected so that it maximizes Rqq* The optimum representation condition 


3R 


00 

00 


(37) 


puts the maximum amount of the reflected radiation in the lowest order mode 
(a different mode and matrix element may be chosen, if desired). It is 
pointed out that (37) does not improve the properties of the reflected 
radiation. It enables one to choose the best representation in terms of 
minimizing the coefficients of the undesired modes for the scattered 
radiation. Once W q is fixed by (37) then the exact location and size of 
the waist(s) for the reflected modes is determined by solving the system of 
equations 


(38a) 


W 


o 


x 2 , 1/2 


o 

The element R^ is identical in both representations since 
order mode Uqq is the same in rectangular and cylindrical coordi 
Performing the integration (29a) yields R^ to first order in £ 


(33b) 


the lowest 
nates. 
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R 00 

R oo = 


\ _ e -(l + a 2 )U 2 ' 


(39) 


Note that the first order term vanishes and the lowest correction is of 
2 

second order in £. . The exact dependence on the mirror size p is 
parametrized by y = p/W Q , while a = W./W q parametrizes the ratio of the 
incoming and scattered radiation spot sizes at the mirror. The optimization 
condition 3 Rqq( 0)/9« = 0 yields, a 2 = 1 + exp[-(l+a 2 )y 2 ][2y 2 a 4 +(2y 2 +l)<x 2 -l]. 
In case that the mirror cross section is much larger than the spot size of 
the incoming mode, y » 1, a -» 1 and the reflected spot size at the mirror 
matches the incoming, W Q = Vj. 

Large mirror size is desired to maximize the total reflection 

coefficient. For incoming radiation of unity electromagnetic flux 
Pj = |c*1 2 = Z l c * p q| 2 = 1» the total reflection coefficient tr^ = P Q /P i 
equals the reflected flux P Q , 

P o * I c °I 2 “ I R ’ c ’| 2 ■ E S '‘pql 2 ’ <«> 

mn pq 


In Fig. 3 we plot for the lowest order incoming mode as a function of 
y' = cos($/2) p/W Q = cos($/2) y. y' parametrizes the size of the mirror 
projection into the plane perpendicular to the incoming radiation 
direction. The incoming radiation has a wavelength X = ly (10~ 4 cm), 

-1 2 

waist Wj = 2x10 cm at distance 1^ = 1.8x10 cm from the mirror and 

radius of curvature (at the mirror) R,, = 8.95xl0^cm. The mirror has a 

radius of curvature R = 8.95xl0^cm, yielding reflected modes of R = 

m o 

3 

8.95x10 (again 1 q and v q depend on the choice of W o ). In Fig. 4 we 
plot the magnitude of the reflection coefficients |Ppq| of the lowest 
order mode (0,0) into the firstr 25 modes (p,q) with p < q < 5, as a 
function of y'. The deflection angle is 90° and the ratio of the spot 
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sizes is 1. Increasing mirror size maximizes the diagonal element and 
minimizes scattering into other modes. The spherical aberration was 
retained inside (23) in evaluating the matrix elements. Its effect is 
small, since for the above parameters £ = 0.28S, and a good agreement 

is observed with the constant A theoretical limit. In particular, the 
dominant off-diagonal terms couple the (0,0) incoming mode to the (1,0), 
(3,0) and (3,2) reflected modes only, according to the selection rules, 
Eqs. A(10). Comparing Figs. 3 and 4 with the next plots shows that the 
relative mirror size to the radiation spot size is the most important 
parameter to determine the reflection into other than the incoming 
modes. 

In Fig. 5 we fix the mirror size y' = 2 and the angle $ = 90° and 

vary the spot size ratio a. The best representation, maximizing R^® 

and minimizing R 00 is obtained at a = 1. However, for small mirror y' = 
pq 

0.66, the maximum for Rqq occurs at a = 0.70 (see Fig. 6). Radiation 
reflected off mirrors smaller than the incoming spot size is best 
described by outgoing modes of reduced spot size V Q < W^. Also note 
from Fig. 6b that for small mirror size the total power reflected into 
the first 25 modes never exceeds 80% of the incoming flux; even with 
many more modes remains less than 1. In Fig. 7 the reflection 
coefficients R^j are plotted as functions of the angle of deflection $ 
for fixed a = 1, y' =2. It is seen that, for sufficiently large 
reflecting surface and good choice of the spot size W q , the reflection 
matrix is not very sensitive to ^ and the off-diagonal terms remain 
small. 

The main conclusions so far are summarized as follows. When the 
mirror size is > 2.5 times the incoming spot size, the fraction of the 

2 

incident power scatterd into different modes is of order £ for 
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k.V^/R m < 1, or (£*)^ for kV^./R ffl > 1. This holds for a wide range of 
deflection angles $. It will be shown in the next section that cross¬ 
polarization effects are of the same order. In most applications both £ 
* -2 

and £ are less than 10 . To this end, scattering losses will be 

smaller than the losses caused by the finite reflectivity (i.e., 
absorption) by the mirror, for most dielectrics. 
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VII. CROSS-POLARIZATION EFFECTS 


The curvature of the mirror surface produces a cross-coupling between 
the transverse components of the incoming and reflected radiation. 
Inserting expressions (21) for the normal unit vector to the mirror inside 
the boundary conditions Eq. (7a), the full source term 

A s = (A X S , A y S , A/) for an incoming wave A 1 = (A^, A y *, 0) is given by 


2 cos 2 r V 


(41) 


i s 

In the above relations, the components of A and A are given in coordinate 
systems aligned with the incoming and outgoing radiation, respectively. 
According to (41) the reflection of a plane polarized wave generates 
components polarized in every direction, including A z . These cross 
polarization effects enter to order £ and result in a small rotation of the 
polarization angle. 

The radiation steming from the A^ s component will propagate 
perpendicularly to the direction of interest z q and escapes the resonator 
as pure reflection loss. The relation between the incoming and reflected 
transverse components, including cross-polarization effects, is now given 
by 



The matrix 3 has been computed in the previous section. Substitution 
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of the additional cross-terms in Eq. (41) inside the propagator integral 
(6) yields 

Q mn = 0 mn = O'™ 

(xy)pq (yx)pq w pq ’ 


fl* 


- i5 mn< x s *y s > 



1 

fl , z ° 2< vv 1 

V.<V> , s >U Da<W 

-O 2 


r,, \LT 



1 >0 2 J 

1 * b. 2 

1 

'V 

[cos H 1 - 4 sin 




m J 


In Gaussian-Hermite representation, we obtain 

W r > • 

wi th 


Q™(1) = Jll C“ e pq tan I G “ 
pq pq 2 pq 


Y o 


pq 


M 


dY H («X) H (Y) Y H (X) H (Y) . 


(44b) 


In Gaussian-Laguerre representation, we have 

, *. 

QPJ(l) = n l sin | J dX Gj q (X) B pq (X), 


where 


i ^ 

2 x ^ 2 x 2 i p r* 2 x^ i.V 


G pm (X) = (a X) X L F (a X) L l, (X) e 


X 


(45a) 


(45b) 
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and 



In both representations, cross polarization effects enter to order £. In 
case of vertical incidence ($ = 0) with arbitrary curvature 1/R m , or plane 
mirror (£ - 1/R m = 0) and arbitrary incidence <J>, Q goes to zero. 

Transverse vector components are reflected independently of each other in 
these two limits. Some of the elements of Q (in both representations) are 
given in Appendix C for large (p » Wj) mirror. 


Acknowledgmen t 


This work was supported by SDIO and managed by SDC. 


28 




Appendix A. Computation of the Hermitian Matrix Elements. 

The integrals (29) will be evaluated here in case the mirror size p is 
much larger than the incoming mode spot size Wj, p cos +/2 » Then the 

limits of the surface integrals can be extended to infinity, and the spot 
size for the outgoing modes V matches that of the incoming at the mirror, 
i.e., a = 1. We use the notation 

+ n = e 2 H n (X), <A1) 

the recurrence relation 

H n '(X) = 2n H nl (X), 
and the orthonormality properties 

I dx *„< x >»„< x > - V»> 


21 t 2 "' 1 " 14 .,™-! * 

(A4) 

' { 2 "‘ Z " !S ., n -2 * 2 "' 1 < 2 "* 1 >"’ S ,,„ * 2 "<"* 2 >' 5 n ,,„, 2 }- < A5 > 

to obtain 

I™ 1 = n 2 m+n n!m! & 8 , (A6) 

pq P* m Q.n 

M™ = - \ <X> f/n 2%!)8 , 

pq Z p, m t ) q, n 


J dX * n (X) X * m (X) a <X> m>n = 

J dX V X > x2 V X > a <x 2 > m,n 


(A2) 


(A3) 
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- T [ 2k<x Vi,. - 2 * <x2> .-iJ( /; 2 ”"'K.» 

- J 1 [2k 2 m m! » p . liB - 2. 2- 1 <.-l)!8 p>m _ 1 ]<2 2 >,,„. <«> 

N “ . 1 * <p*q*DjV[] <* 8 > 

Inserting expressions (A3)-(A5) into Eqs. (A6)-(A8) we obtain 

„mn , nx c i (A9) 

R p, <0) ■ V* S q.n’ 

«-»> - 


(A10) 


«»f * 7 r ( 5 r) 3/2 - W ♦ W /I ls.-1 

- g m i/2 - * k pfl 1/l & * 

- L piliap^) 1 ' 2 W3 . L P^ > ) 1/ V»-3} 

— [««*!,;$. (*,.!>&] [&T^ * 

/2 L i i 0 ° 
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R is diagonal to zeroth order. The lowest order correction couples m with 

m±l, m+3 in the X-direction and n with n, n±2 in the Y-direction. The 

reflection matrix is not symmetric, R™ 1 * • Also, it is not invariant 

pq mn 

to interchanging X and Y. This means that the modes u mn ( x >y) and u nm ( x >y) 
with m * n are reflected differently. 
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Appendix B. Computation of the First Matrix Elements in Laguerre 


Representation. 

Representation using Gaussian-Laguerre modes may be advantageous in 
numerical simulations because fever Laguerre modes than Hermite modes are 
required to represent close-to-axisymmetric radiation profiles with the 
same accuracy. However, the computation of Eqs. (33a) to (33e) is not so 
straightforward. The integrations (33a) and (33b) for I ptl and K pq over the 
polar angle 0 s involve trigonometric functions of complicated arguments 
e^Gg) and 6 0 (© s )» given implicitly by Eq. (31b). There is no simple 
recurrence formula for this calculation. The first few elements are 
computed here by expansions in powers of r s /R m < C. - £. Substituting from 
(31b) inside (33) and renormalizing = X/( l-sin 2 4>/2sin 2 e), one 

obtains, to first order in £., 


(BO) 
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D 1-1 (X) = 0. 


(B5) 


We only need E pq to zeroth order in £, obtaining 


U (X) = 


l m -1 


t ,in 2 e 


cos ? sin 


^1-sin 2 | sin 2 9j 


and 

E pq = 0 + 0(£.) for (p,q) t (-1,0), (0-1). (B7) 

The integrals (Bl) - (B7) are evaluated using the formula 



o 


u v 
sin x cos x 

(l-k 2 sin 2 x) p 


¥) p ( p * * ± t 1 ' v2 \ 


(B8) 


where B(p,q) = r(p)r(q)/T(p+q), T is the factorial function and F is the 
hypergeometric function. The radial integrations for U, V and W are 
performed directly, using the expressions L p (x) for the Laguerre functions 
and the identities 


J e“ x x" 1/2 dx = /it, 
o 

J e -x x n dx = nl, 
o 

J e _x x n+1/2 dx = | • | ...(n+l/2)/n. (B9) 

o 

Again, we extend the limits of integration to infinity assuming p cos($/2) 
» and a = 1. The zeroth order contribution is given by 

O o > = 

rJJ(O) = 1, 
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R 1 ~ 1 (0) = cos3 2 F ( 2 .|. 2 » sin2 f )’ 

r}}( 0) = cos | F(2,|,2,sin 2 f), 
rJJ< 0) = cos | F(2,2,sin 2 f), 

R q q(O) = cos 2 | F(2,|,2,sin 2 f). 

The first order corrections in f. are given by 

R 00 (1 > • - R 00 (1 > * (I ) 1 2 I sl " 2 F ( 2 ’J' 2 ' 511,2 !)’ 
R 0i a > ■ - R ?J (1 > ■ (f) 1/2 t6 sl " 2 F ( 2 ’?' 2 - 5i " 2 £)• 

R “(D . R°°<1) - 0, 

r“< 1) - -R^U) - (f) ^sin * F(2,i 2, sin 2 *), 


(BIO) 


R ll (1) = ~ R 11 (1) = - (f) ff sin 2 F ( 2 ’?’ 2 ’ si ” 2 £)’ (B11) 


r -J 2<1) = sin | |± |(f) | F (2,i 2, sin 2 f) - cos 2 f F ( 2 ,2, sin 2 

1 1 c„s 2 ^g,|,2, si„ 2 

{* [d) 1/2 16 F ( 2 >|' 2 ’ 

g-,|,2, si„ 2 |)[2 - i 2- 122 ]}. 


R'JJd) « sin | 


cos 2 | F f2,^,2, sin 2 


,s 2 f FfU.S 
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R ;"(!) = sin ; 


± \ COS 2 | F 


| t [g) 1/2 gp( 2 ,i 2 , sin 2 


I 2 £)- cos 2 | P(2,|,2, sin 2 |) 


The (-) sign and the lowest row inside the last square bracket in (B12) 
correspond to exchanging indices, 


of order £. or higher for every m, n, 
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Appendix C. Computation of the Cross-Polarization Matrix Elements 

We compute here some of the first order cross-polarization matrix 
elements in case of large mirror size p » and a ■* 1. In the Gaussian- 
Hermite representation ve find from (44b), using the notation of Appendix 
A, 

G®" = 2 m m! /n <Y> 6 

pq q,n Pim 

= n 2 m m! J2 n_1 n! S , + 2 n (n+1)! $ ,1s (Cl) 

1 q,n-l q,n+l) p,m, ' ' 

yielding 

0“ (1) - t .in I {/I 5 P ,„. (C2) 

In Gaussian-Laguerre representation we only have to compute B pc| (X), 
Eq. (45c), to zeroth order in £. Applying the methods of Appendix B, ve 
find 

B pq = 0(0, if p,q * (1,0), (0,1), 


o (l-«in 2 f S 


(C3) 


Noting that G mn (X) is the same as w j^( x )> Eq. (33e), and inserting (C3) and 
(33e) inside (45a), we obtain 


Q\° 0 (1) = Q° 0 \ (1) = 0, 

Q oo = Q oo = 1 sin I F d’ b 2> sin2 I)’ 

Q 01 (1) = Q 10 (1) = sin 2 p (|» b 2> si " 2 !)’ 

Q 11 (1) = Q 11 (1) " /2 1 sin 2 F d» b 2 » sin2 !)• <C4) 
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Appendix D. Small Aperture Limit 

We have seen in Sec. V that in case of mirror surface S^ with an 
aperture of area S2 the reflection matrix is given by 

R(S;V o ) = R(S 1? W o ) - R(S 2 ;W o ) (Dl) 

In case that S 2 « the spot size W q optimizing the representation 
for the scattered radiation will be determined predominantly by the surface 
Sj. Thus, the formula (23) with W q given from 

3R(S.;a) 

Bee - = 0 . < D2 > 

can be used for the modal decomposition of the scattered radiation. 
According to Eqs. (28a) and (29a) for the Hermite representation, and Eqs. 
(32a) and (33c) for the Laguerre representation, the lowest order 

2 

contribution from a small aperture p 2 « W Q scales as R(S 2 ;W o ) ~ £ . 

In some cases, however, it is important to know the total radiation 
diffracted through a small hole, rather than the modal decomposition. In 
case of small apperture p 2 » 

p 2 2 « l Q k -1 or X » Pj/Iq. (D3) 

where 1 q - z is the distance of the observation point from the mirror, the 
paraxial approximation, Eq. (8b) is taken one step further, setting 




Substituting (D4) inside (8a) we obtain the "far field" limit of the 
diffracted radiation 


A°(x,y, z=0) = ^-Jdx o Jdy o A 1 (x o ,y o )(n-8r) 

“P { - ^ [(* 2+ v 2 ) - 2 «„ - 2 W„]}’ < D5 > 

also known as Fraunhoffer diffraction. The condition (D3) can only be 

valid for appertures much smaller than the spot size W q at the mirror, p 2 

« V Q , since for p 2 - W Q (D3) is violated, kp 2 2 - kV Q 2 ~ kw Q 2 (l+l o 2 /b o 2 ) 

~ b (1+1 2 /b 2 ) > 1 . Neglecting terms of order kx 2 /l - kp 0 2 /l << 1 
o oo o mo2o 

means that terms of order x /I , x /I « 1/kx , where kx >1, must also 
s m s o s s 

be neglected. The source term can be written as A*[x.,y^] = A*[x g ,y s ]. 
Rescaling variables to 

K x . cos K y . ^ (D6) 

we obtain 
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2 2* 2 
x s cos 2 + y s 


cm I KK "p[ /? 57 cos fhl* 57] 


i(m-n)6 .m+n 2 


The zeroth order contribution in x /V. « 1 yields 


ika ik 1 + 

A°(x,y,0) = o' e 1 0 


X cos | Jdx g Jdy s J o (Kx g ) J o (Ky g )e 


x s 2 cos 2 I . y s 2 


The waist size w f for the Fraunhoffer inodes is given by the zeros of the 
Bessel functions 


K(v f ) x - — p» ~ 2n. 


Therefore, the diffraction angle 9 f is 


e f ■ A - K 


The requirement 9^ « 1 for the validity of the paraxial approximation puts 
a lower limit in the aperture size Pj 
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In case the aperture size is of the order of the vavelength X the 
scattered vavefunctions are spherical rather than Gaussian. Because the 
overall effect of a scatterer with size P 2 ~ X is very small, the familiar 
from quantum mechanics Born approximation, involving perturbation theory, 
is applicable in that case. 
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Figure 2 Reflection geometry. 
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3 Plot of the total reflection coefficient for the lowest o 
mode as a function of the mirror size u' for $ = 90°. The 
radiation has wavelength A = 10'^cm, waist w i = 2xl0 -1 cm at 
2 

distance 1. = 1.8x10 cm from the mirror and radius of curvati 


8.95x10 cm. 
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7 Plots of the reflection matrix elements 


against the angle of 


deflection $ for y' = 2 and a = 1. 
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